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Abstract - To consider relation between percolation criticality and network-growth rule in an 
evolving network model, we propose a general network-growth rule which includes the conventional 
network-growth rule and the rule proposed by Achlioptas et al. very recently [Science 323 (2009) 
1453]. We introduce the generalized parameter q which characterizes how to evolve a network. 
We obtain size dependency of percolation step and fractal dimension of percolated cluster at the 
percolation step for several g's. Strong size dependency of percolation step and sudden change of 
fractal dimension are observed for intermediate q. From these facts, we find that the percolation 
criticality changes between at the conventional network-growth rule and at the rule proposed by 
Achlioptas et al. 
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■ Percolation transition has been a hot topic not only 
In statistical physics but also in a wide area of science 
[T}{3]. In materials science, theory of percolation tran- 
sition is used to investigate relation between connectivity 
and properties of materials such as electric conductivity in 
Alloys, magnetic phase transition in diluted ferromagnets, 
and sol-gel transition in polymers. In natural science, pcr- 
tolation transition relates to dynamical features in nature 
such as spreading wildfire and spreading epidemic. Al- 
though these phenomena arc non-linear non-equilibrium 
behaviors, we can consider these effects by just focusing 
(bn connectivity of composing elements in the network. In 
information science and engineering, theory of percolation 
transition can also adopt for dynamic nature of evolving 
network systems such as complex network. A curious ex- 
ample is to investigate relationship between percolation 
transition and PageRank which is an algorithm for decid- 
ing importance of website and is used in search engine 
[3]. Applicable scope of percolation theory has been ex- 
tended day by day. On the other hand, it seems that 
nature of percolation transition such as relationship be- 
tween the criticality and the spatial dimension has been 
well-established in terms of statistical physics. 

In 2009, Achlioptas et al. proposed a novel type of 



network-growth rule and studied nature of percolation 
transition [S]. In the conventional network-growth rule, we 
randomly select a bond between elements which do not be- 
long to the same cluster and connect the bond definitely. 
On the other hand, in the rule proposed by Achlioptas 
et al., wc randomly select two bonds between elements 
where the elements adjoining each bond do not belong to 
the same cluster. After that we connect the bond where 
the product of sizes of two clusters which arc contacted by 
the bond is smaller than the other. Wc refer to this rule as 
Achlioptas product rule. They calculated time evolution 
of the maximum size of clusters on system without spa- 
tial structure. Nature of time evolution resulting from the 
Achlioptas product rule is completely different from that 
obtained by the conventional network-growth rule. In fact, 
Achlioptas et al. claimed that a discontinuous phase tran- 
sition takes place in network-growth model adopted the 
Achlioptas product rule, whereas a network-growth model 
using conventional rule exhibits a continuous phase tran- 
sition. Since the size of the maximum clusters increases 
explosively against time in the Achlioptas product rule, 
the percolation transition is called explosive percolation 
transition. Some researchers have confirmed occurrence 
of explosive percolation transition [HIIZIj whereas a num- 
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Fig. 1: (Color online) Graphical representation of our pro- 
posed rule on the square lattice. The numbers in the diamonds 
represent the state of vertices. The transparent gray diamonds 
denote the isolated clusters, (a) The bold lines denote ran- 
domly chosen edges dj and eu- In this case, the sums of the 
clusters are n(l) -I- n(3) = 12 and n(4) -I- n(6) — 13. We select 
eij{eki) as connecting edge with the probability Wij{wki)- In 
the random rule {q = 0), we select e^j or e^i with the same 
probability 1/2. In the Achlioptas sum rule {q = +oo), we def- 
initely select Eij whereas we select eki in the inverse Achlioptas 
sum rule {q — — oo). (b) Configuration after the edge e^j is 
selected and connected, (c) Configuration after the edge e^; is 
selected and connected. 



bcr of researchers have insisted that the explosive per- 
colation transition is actually continuous [5Hl2|. After 
this breakthrough study by Achlioptas et al., a number 
of researchers have thoroughly investigated properties of 
explosive percolation transition on finite-dimensional lat- 
tices (regular lattices), scale- free networks, and random 
graphs by adopting the Achlioptas product rule [T3HIS]- 
Moreover, a couple of extended network-growth rules have 
been proposed and nature of percolation transitions using 
these rules have been studied [T7H20] . The aim of these 
studies is to propose a rule which makes the percolation 
transition more explosively. In contrast, our purpose of 
this letter is to investigate relation between percolation 
criticality and network-growth rule. To study this prob- 
lem in a clear manner, we adopt the sum rule as network- 
growth rule in this study. In the sum rule, we randomly 
select two bonds between elements which do not belong to 
the same cluster and connect a bond by comparing sums 
of the size of clusters whereas in most preceding studies, 
products of the size of clusters are compared. When we use 



the sum rule, the number of elements obeys additivity be- 
tween before connection of elements and after connection 
of elements because of the definition. In other words, the 
small cluster is definitely constructed in network-growth 
process based on the Achlioptas sum rule in which we 
connect the bond where the sum of sizes of two clusters 
is smaller than the other. Then we can figure out what 
happens clearly when we change how to evolve a network. 
The sum rule is the simplest and most appropriate rule 
to control the percolation phenomena such as the change 
of percolation step depending on the network-growth rule. 
To control a network-growth rule, we introduce a new pa- 
rameter q which characterizes how to evolve a network. It 
should be noted that this is the first study where the sum 
rule is adopted for two-dimensional lattices. Our proposed 
rule includes the conventional rule and the Achlioptas sum 
rule as shown below. 

Before we explain our proposed rule, we define network- 
growth model and notations in a general way. Let V and 
-E be a set of vertices (points) and edges (sides), respec- 
tively. The i-th vertex is denoted by Vi € V, and € E 
represents the edge between Vi and Vj . The state of the 
i-th vertex is represented by Oi [oi = {1, • • • , N}), where 
N denotes the number of vertices. The state of the edge 
between the i-th and the j-th vertices is expressed by Tij 
{Tij = {0, 1}). We assume that all edges are non-directed 
i.e. Tij ~ Tji for Vi, J. When the edge e^j is not connected, 
Tij =0, whereas when the edge is connected, Tij = 1. 
Cluster is defined as a set of vertices which have the same 
state. The number of elements in the cluster where Vi 
belongs to is expressed by n{ai). 

Here we explain our proposed network-growth rule. We 
show that our rule can describe the conventional network- 
growth rule which is referred to as "random rule" and 
the Achlioptas sum rule. The procedure of our proposed 
network-growth rule is as follows: 

Step 1 The initial state is set to be ct; = i (i = 1, ■ • • , N) 
for Vi, in other words, all elements belong to different 
clusters, i.e. n(ai) = 1 for Vi and Tij = for Vi, j. 

Step 2 We randomly choose two different edges e,;j ^ Cki 
satisfying conditions such that = Tki =0, ai ^ Uj 
and (Tfc 7^ cr;. 

Step 3 We connect e,;j with the probability Wij defined 

by 

g-9[ri((Ti)+Tl((Tj)] 

•= c-9[n(^.)+"(<^.)] +e-9["K)+"('^')] ' '^^^ 

whereas we connect Cfe/ with the probability w-ki 
1 — Wij. After we connect eij{eki), the states of one 
of the clusters where Vj{vi) belongs to are changed 
to Uiiuk). When we connect an edge, time advances 
from T to T+ 1. 

Step 4 We repeat step 2 and step 3 until all of the ele- 
ments belong to the same cluster. 
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Here we assume that all clusters are never separated. In 
this network-growth rule, the number of clusters decreases 
one by one in each time. In fact, the number of clusters 
at time risiV-rforO<T<7V-l. The graphical 
representation of the above procedure on the square lattice 
is summarized in Fig. [TJ 

The introduced parameter q in Eq. ([1]) is a generalized 
parameter and characterizes the network-growth rule. Our 
proposed rule for g = is equivalent to the random rule 
as follows. In the random rule we randomly choose an 
edge Cij such that cr^ ^ aj and = and connect the 
edge Bij. In Eq. ^ ioi q ~ 0, the probabilities are the 
same Wij = Wki = 1/2, and thus it is equivalent to select 
the connecting edge randomly. On the other hand, our 
rule for q = +oo realizes the Achlioptas sum rule. In this 
rule, we randomly choose two edges and Cki such that 
the conditions which stated in step 2 arc satisfied. We 
compare the sums n{ai) + n(crj) and n{ak) + n{<Ji) and 
connect the edge where sum is smaller than the other. 
Actually, for q = +cx) in Eq. (IT|), wtj = 1 and wu = when 
n((Ti) + 7i(fTj) < n[ak) + n{(7i) and vice versa. Hereafter we 
refer to our rule for q = — oo as "inverse Achlioptas sum 
rule". In this case we select the connecting edge where 
the sum of elements is larger than the other. In this way, 
network-growth rule can be controlled by the generalized 
parameter q. Thus, we can explore relationship between 
network-growth rule and nature of percolation transition 
through the generalized parameter q. Notice that because 
of the conditions imposed in step 2, a cluster has no loop as 
well as "loopless" percolation [21]. This treatment is the 
same as the Achlioptas product rule process considered by 
Ziff [TS]. 

To investigate the network-growth dependencies of na- 
ture of percolation transition, its criticality, and fractal 
dimension of the percolated cluster, we study the network- 
growth model on Lx L{= N) two-dimensional square lat- 
tice with open boundary condition. The coordinate of 
Vi is represented by the position vector r; := (xi,yi) for 
1 < Xi,yi < L. 

In order to understand g-dependcncy of network evolu- 
tion quantitatively, we calculate the dynamics of the num- 
ber of elements in the maximum cluster rimax defined by 
'^max := max{n(Q;)|l <a< N}. Here n{oi) represents the 
number of elements such that ct; = a (1 < i < TV). The 
value of rimax characterizes the connectivity of network 
and is often calculated in studies of conventional percola- 
tion. As stated above, since we assume that clusters are 
never divided, rimax monotonically increases with time. 
Figure [5] shows time development of the density of ele- 
ments in the maximum cluster nmax/-^ for q = —oo (in- 
verse Achlioptas sum rule), —1, — 10~^, —10^^, —10^'^, 
-10"'', (random rule), 2.5 x 10-^ 5 x 10-^ 10^"^, 
2x10""*, 10~^, 10"^, 10"\ and +oo (Achlioptas sum rule) 
on 256 x 256 square lattice. These results are obtained by 
averaging out the results calculated over 1024 samples. 
Here we define t as the normalized time t := T/N. The 
step of rise monotonically increases as q increases as shown 
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Fig. 2: (Color online) Time evolution of the density of ele- 
ments in the maximum cluster nn-iax/N on 256 x 256 square 
lattice for q = —oo (inverse Achlioptas sum rule), —1, — 10~^, 
-10"^ -10"^ -10"", (random rule), 2.5 x 10"^ 5 x 10"^ 
10"*, 2 X 10"*, 10"^, 10"^, 10"\ and +oo (Achlioptas sum 
rule) from left to right. These results are obtained by averaging 
out 1024 independent samples. Since the error bars are smaller 
than the widths of lines, they are omitted. 



in Fig. [5] Thus, we can control network-growth speed 
by tuning network-growth rule. In contrast, the gradient 
at the step when n^^x/N reaches near 1/2 behaves non- 
monotonically against the parameter q. This fact relates 
to the size dependency of the percolation step as explained 
later. 

So far, we constructed a general network-growth rule by 
introducing the generalized parameter q and showed the 
g-dependency of dynamical nature of network evolution. 
Next we concentrate on the (j-dependency of percolation 
point. In this study we define percolation as follows. If 
there is a cluster that spreads from left-side to right-side 
or from top to bottom, we call it percolated cluster. This 
is a typical definition in the percolation theory. Here we 
consider size dependency of the percolation step tp{L) for 
several g's. The percolation step is defined by the step 
when a cluster becomes first percolated. The value of 
tp{L) is obtained by averaging out the obtained individual 
percolation step from 1024 samples. Figure |3] (a) shows 
size dependency of percolation step tp{L) for several g's 
which are the same values using Fig. [5J As q increases, 
the percolation step increases monotonically. This fact 
is consistent with the result shown in Fig. [21 The per- 
colation step tp{L) strongly depends on the lattice size 
around q ~ 10~^, which is the reason why the gradient 
of fimax/A^ at the step when rimax/A^ reaches near 1/2 be- 
haves non-monotonic against q in Fig. [21 Figure [31 (b) 
shows g-dependency of the percolation step for L ~ 256. 
At q ~ 10"^, this value suddenly increases as q increases. 
Next we study g-dependency of geometric aspect of the 
percolated cluster at the percolation step. To consider it, 
we calculate the ratio between the number of elements in 
the percolated cluster rip and that of surface elements in 
the percolated cluster ng. The surface element means the 
element that adjoins the other clusters or an isolated ele- 
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the fractal dimension V is adopted as follows: 
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Fig. 3: (Color online) (a) Size dependency of the percola- 
tion step tp{L) for q — — oo (inverse Aehlioptas sum rule), —1, 
-10~\ -10"^ -10"^ -10""*, (random rule), 2.5 x 10"^ 
5 X 10"^ 10"*, 2 X 10"*, 10"^ 10~^ 10"\ and +oo (Aehliop- 
tas sum rule) from bottom to top. (b) q-dependency of the 
percolation step tp{L) for L = 256. (c) q-dependency of the 
ratio ris/rip for L = 256. These results are obtained by aver- 
aging 1024 independent samples. Since the error bars of these 
results are smaller than the symbol sizes, they are omitted. 



ment. Figure[3](c) shows dependency of the ratio ns/n-p 
for L = 256. At q c± 10^^, this value suddenly decreases 
as q increases. In other words, geometry of the perco- 
lated cluster abruptly changes at near q ~ 10"^. These 
obtained results indicate that the percolation criticality 
around q ~ 10"^ is different from that of other regions 
including the random rule and the Aehlioptas sum rule. 

In order to investigate relation between the geome- 
try of percolated cluster at the percolation step and the 
network-growth rule more quantitatively, we consider q- 
dependency of fractal dimension. The fractal dimension 
T) can be calculated by the radius gyration TZ and the num- 
ber of elements in the percolated cluster n-p. The radius 
gyration of the a-th cluster 'R.{a) is defined as 



n{a) := 



^ |r, -ro(a)|2. 



n{a) 



i s.t.cTi —a 



ro(a) 



E 



(2) 



(3) 



where the summation takes over vertices such that ai = a 
and ro(a) represents the position vector of the gravity 
point of the a-th cluster. In this study the definition of 
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where TZp denotes the radius gyration of the percolated 
cluster at the percolation step. The fractal dimension 
quantitatively characterizes geometric properties of frac- 
tal systems and is often used in analysis of the fractal 
geometries. The fractal dimension displays a criticality in 
percolation transition. When D ~ d, where d is a spatial 
dimension, the geometry has no fractality. If < 23 < d, 
on the other hand, the geometry of the percolated cluster 
has fractality as well as conventional self-similar structure. 
The upper panels of Fig. H] show snapshots of the perco- 
lated cluster and those of the second-largest cluster for 
q = — oo (inverse Aehlioptas sum rule), —10^^, (random 
rule), 10"^, 10"^, and, -|-oo (Aehlioptas sum rule) from left 
to right. The corresponding radius gyration dependency of 
the number of elements in the percolated cluster is shown 
in the lower panels of Fig. 21 which are obtained by cal- 
culation from L = 64 to L = 1280. The dotted lines are 
obtained by least-squares estimation using Eq. (jj]). Since 
these graphs are double logarithmic plots, the gradients 
of these fitting curves express the fractal dimension. The 
values of fractal dimensions are denoted in Fig.jH The ob- 
tained fractal dimensions for all negative q are the same 
as that for q = (random rule) T) ~ 1.88 within the error 
bar. Moreover, for small q{< 10"^), the fractal dimension 
is also the same as that for g = within the error bar. It 
should be noted that the fractal dimension for this region 
is the same value as the conventional percolation transi- 
tion on square lattice [3]. On the other hand, the fractal 
dimensions for positive value g(> 10"^) are the same as 
that for the Aehlioptas sum rule V ~ 1.97 within the error 
bar. The fractal dimensions of the random and Aehlioptas 
sum rule are clearly distinguishable. Here we should no- 
tice that the fractal dimension for Aehlioptas product rule 
is the same value as in the case of Aehlioptas sum rule. 
Then qualitative similar behaviors are expected even when 
we adopt the rule based on product rule as most preceding 
studies. Thus we conclude there are at least two regions 
of q in terms of fractal dimension, and between at these 
regions the percolation criticality changes. It is an open 
problem to investigate the change of the fractal dimension 
in the intermediate q in detail. 

In this study we proposed a novel network-growth rule 
in which the introduced parameter q assigns the rule. 
Since our rule includes the Aehlioptas rule which was in- 
troduced in 2009 15] and the conventional network-growth 
rule, our rule can be regarded as a general network-growth 
rule. We studied time evolution of the number of maxi- 
mum cluster. We concentrated on the rule dependency 
of percolation criticality which is characterized by size- 
dependency of the percolation step and the fractal dimen- 
sion. Strong size-dependency of percolation step in the 
case of g ~ 10~* were observed. Fractal dimensions for 
several q were also calculated. The fractal dimensions at 
all negative q and the small positive value (g < 10^^) 
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Fig. 4: (Color online) (Upper panels) Snapshots at the percolation step for q = — oo (inverse Achlioptas sum rule), -10-^ 
(random rule), 10~^, 10~^, and, +cxd (Achlioptas sum rule). The dark and light colors depict the percolated cluster and the 
second-largest cluster, respectively. (Lower panels) Number of elements in the percolated cluster rip as a function of radius 
gyration TZp for corresponding q. The dotted lines are obtained by least-squares estimation using Eq. (|4]| and the fractal 
dimensions are denoted. 



are V ~ 1.88. On the other hand, the fractal dimen- 
sions for 10^^ < 9 < +00 are V ~ 1.97. Thus, we con- 
clude that the criticality of percolation transition changes 
around q 10~^ by considering network-growth rule in 
full detail. 

In this study we focused on the ease for two-dimensional 
square lattice. To investigate the relation between the spa- 
tial dimension and the percolation criticality for our pro- 
posed rule is a remaining problem. Moreover, our rule is 
a general rule for many network-growth problem and en- 
ables us to design percolation criticality. Then we strongly 
believe that our rule will provide a greater understanding 
of percolation transition and will be applied for network- 
growth phenomena in nature and information technology. 
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